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1. Introduction 

Analysis of Deep-inelastic scattering (DIS) cross sections is generally performed through 
global fits of scale-dependent quark and gluon distribution functions fi{x,Q^), i = q,q,g. The 
underlying theoretical framework of such analysis is based on the collinear factorization theorem 
[1], which organizes the computation of DIS structure functions F 2 ^l{x,Q^) into the convolution 
of short-distance Wilson coefficients and long-distance parton distribution functions (PDFs) (see 
e.g. [2] for a recent review). In order to formulate such a theorem it is necessary to introduce a 
factorization scale Hf which separates long- and short-distance physics. Independence of physical 
observables on /i/ allows then to derive renormalization group equations (RGBs) which govern 
the scale dependence of PDFs, known as the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equations. Since there is an infinite number of different ways to realize factorization, one is left with 
an additional choice of the factorization scheme with the MS prescription the generally adapted 
one. For observables, such as DIS structure functions, any residual dependence on factorization 
scheme and - scale jlf is suppressed by an additional power of a^, i.e., is formally one order higher 
in the perturbative expansion but not necessarily numerically small. 

As an alternative to this conventional treatment it is also possible to formulate QCD scale 
evolution equations directly for observables without referTing to auxiliary, convention-dependent 
PDFs, which avoids introduction of an artificial facforizalion scheme and -scale dependence al¬ 
together [3]. The framework is suited best for theoretical analyses of DIS data; in particular one 
remains in this case with the renormalization scale as the only theoretical ambiguity. Since the¬ 
ory uncertainties are in this way reduced to a minimum, it is therefore this ‘physical’ formulation 
of DGLAP evolution which is most suitable for extractions of from inclusive DIS data (for a 
first study see e.g. [4]). Moreover, physical evolution allows for the most stringent tests of DGLAP 
evolution itself. This is of particular interest for regions of phase space where a breakdown of 
collinear factorization is expected, such as the limit x —0 of DIS structure functions, where large 
parton densities eventually saturate and require a more complete description including terms usu¬ 
ally suppressed by powers of Q^. A potential application of physical evolution for such studies is 
illustrated and discussed in Fig. 1. The outline of these proceedings is as follows: in Sec. 2 we 
provide some details on the derivation and definition of physical evolution kernels while Sec. 3 
presents an analysis of the remaining renormalization scale dependence of physical evolution up to 
next-to-next-to-leading order (NNLO) in for large input scales = 30GeV^. For more details 
we refer the interested reader to [5]. 

2. Physical evolution kernel 

To define physical evolution kernels we use that the x-space convolutions of coefficient func¬ 
tions and PDFs turn into products in conjugate Mellin space, a{n) = dxx"^^a(x). Moments of 

DIS structure functions F/(x, Q^) can then be expressed as 



(2.1) 
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Figure 1: Gold structure functions F 2 and Fi have been calculated at ^ = 2GeV^ from the bCGC-model, 
using the fit [6] to HERA proton DIS data, combined with a scaling ^ The result has been fitted 

and used as input for physical DGLAP evolution. The plots show a comparison of physical DGLAP evolution 
at leading (LO) and next-to-leading (NLO) order for the doublet [F 2 ,Ff) from Qq = 2GeV^ —>■ 2^ = lOGeV^ 
and the corresponding bCGC result at = lOGeV^. 


The sum runs over all contributing quark flavors and the gluon, each represented by a PDF fj^. The 
non-perturbative PDFs obey the DGLAP evolution equations 


dfk{n,lX^) 

dXnfL^ 


l " 


( 2 . 2 ) 


while coefficient functions C]^k [2, 8, 9] and splitting kernels P^i [10, 11, 12] can be calculated in 
perturbative QCD and exhibit the following expansion in 


m=0 
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^bk=L[jz) Q/(n), (2.3) 
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where tuq depends on the first non-vanishing order in Us in the expansion for the observable under 
consideration, e.g. mo = 0 for F 2 and mo = 1 for The DGLAP evolution equation are formulated 
as u/ — 1 evolution equations for the different non-singlet quark flavor combinations and a 2 x 2 
matrix valued evolution equation, which evolves the flavor singlet vector (r,g); the gluon 

distribution and r(n,/r^) = + the quark flavor singlet. In the following 

we concentrate ourselves on the flavor singlet sector only; for the physical evolution of the non¬ 
singlet sector see e.g. [4]. Using that any doublet of flavor singlet observables F = {Fa,Fb) is 
related to the flavor singlet vector (r,g) through a coefficient matrix C 


C = 


( CAq ^Ag 
Csq Cjjg 


P = I 

\ ^gg 


(2.4) 


and introducing further a corresponding 2x2 matrix P for the matrix-valued kernel of the DGLAP 
evolution in the flavor singlet sector, one finds in a straight forward manner 


dF{n,Q^) 

dlnQ'^ 



—c-i 

das 


+ CPC 


■F. 


(2.5) 


where we made in addition use of the RG equation of governed by the QCD beta function 


dasia) , n, ^ 

The resulting physical evolution kernels 
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( 2 . 6 ) 


, (2.7) 


are independent of factorization scheme and -scale [3] with the renormalization scale Pr as their 
only remaining scale ambiguity at finite perturbative order. 


3. Renormalization scale dependence 

To determine the renormalization scale dependence of physical evolution kernels we use ker¬ 
nels calculated at Pr = Q, see e.g. [5], and recover their full renormalization scale dependence 
using the same prescription as for conventional DGLAP splitting kernels i.e. by Taylor expanding 
cCs(Q^) in terms of as{p})\ see e.g. [13] for a discussion in the case of splitting kernels. For the 
following numerical study we fix 71/- = 4 and use a realisfic toy input at = 30GeV^ for quark 
singlet and gluon distribution (taken from [14]), 

xL{x) = 0.6x-‘’- 3(1 -x)3-5(l +5x°'^) 

xg{x) = L6x^'’'^(l —x)"^'^(l — 0.6x°'^) , (3.1) 

from which we calculate structure functions using LO coefficients, independent of the actual stud¬ 
ied perturbative order. The strong coupling is obtained from solving Eq. (2.6) with the QCD beta 
function truncated at the corresponding perturbative order and initial conditions fixed to as{Qo) = 
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Figure 2: Renormalization scale dependence of the doublet {F 2 ^Fl) 
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Figure 3: Renormalization scale dependence of the doublet {F 2 ,Fs) 


0.2. We study as examples the flavor singlet sector of the doublets {F 2 ,Fl) and {F 2 ,Fs) with Fs the 
F 2 scaling violations 


Fs{x,Q^) 


dF2{x,Q^) 

\nQ^ 


(3.2) 


Our result are depicted in Fig. 2 and Fig. 3 and reveal a very small dependence on the chosen 
renormalization scale at NNLO, with the doublet (^ 2 ,^ 5 ) slightly less sensitive to the variation of 
the renormalization scale. Note that the variations have been performed over a very wide range i.e. 
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/x^o/GeV^ G [6,150] and /i^/GeV^ G [20,500] respectively. 

In conclusion we find a very mild dependence on the renormalization scale, if the initial scale for 
DGLAP evolution is rather large. For the case of small initial scales we refer to the discussion in 
[5]. 
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